Mon. Not. R. Astron. Soc. 000, ITHHl (2011) Printed 25 January 2012 (MN MfeK style file v2.2) 



EVN observations of the farthest and brightest ULIRGs in 
the local Universe: the case of IRAS 23365+3604 



(n: 

o- 

<n: 

a- 
»— »; 

(N 



Cristina Romero-Canizales, 1 - 2 * Miguel Angel Perez-Torres, 1 & Antxon Alberdi 1 

1 Instituto de Astrofisica de Andalucia - CSIC, PO Box 3004, 18080 Granada, Spain 

2 Tuorla Observatory, Department of Physics and Astronomy, University of Turku, Vdisdldntie 20, FI-21500 Piikkid, Finland 



Accepted 2012 January 24. Received 2012 January 24; in original form 2011 November 4 



o 

u 

Of 

6 



O 

in 
o 

(N 



ABSTRACT 

We present high-resolution, high-sensitivity radio images of the ultra- luminous infrared 
galaxy (ULIRG) IRAS 23365+3604. We performed contemporaneous observations at 
1.7 and 5.0 GHz, in three epochs separated by one year from each other, with the 
European very long baseline interferometry Network (EVN). We also present comple- 
mentary Multi-Element Radio Linked Interferometry Network (MERLIN) at 1.6 and 
5.0 GHz, and archival Very Large Array (VLA) data, taken at 1.4 and 4.9 GHz. We 
find that the emission at ^5.0 GHz remains quite compact as seen at different resolu- 
tions, whereas at ^1.7 GHz, high resolution imaging reveals some extended structure. 
The nuclear region has an approximate linear size of 200 pc and shows the presence 
of two main emission components: i) one with a composite spectrum due to ongoing 
non-thermal activity (probably due to recently exploded supernovae and AGN activ- 
ity), ii) another one with a steep spectrum, likely dominated by an old population of 
radio emitters, such as supernova remnants (SNRs). Radiative losses are important, so 
re-acceleration or replenishment of new electrons is necessary. We estimate a magnetic 
field strength of ~ 18 /j,G at galactic, and ~ 175 fiG at nuclear scales, which are typical 
for galaxies in advanced mergers. 
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1 INTRODUCTION 

Galaxies with very high infrared (IR) luminosities (Lir = 
L[8-1000 /iin] > 10 12 Lq) known as ultra luminous IR 
galaxies (ULIRGs), dominate the IR backgr ound and the 
star formation rate (SFR) density at z~2 |Caputi et al.l 
120071 ). Although uncommon at lower redshifts, the pres- 
ence of ULIRGs in the local Universe offers the opportunity 
of investigating their parsec scale structure, while profiting 
from the high angular resolution provided by current in- 
strumentation. Studying ULIRGs in the local Universe is of 
great importance since it can aid to understand their high- 
redshift analogu es which dominate the sub-mm sky (see e.g., 
iLillv et al.l[l999l h 

It is thought that ULIRGs represent a key stage in 
the formation of optical qua si-stellar objects (Q SOs) and 
powerful radio galaxies (e.g.. ISanders et al1ll988l ). A study 
based on HST observations and TV-body simulations point 
to diverse e volutionary paths, not necessarily including a 
QSO phase ()Farrah et al.l [200ll k There is however a gen- 
eral agreement on gas-rich galaxy merging as the origin of 
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ULIRGs l|Sanders fc Mirabellll996l h and on the ubiquity of 
enhanced star-formation, which can be found in combina- 
tion with diff erent flavours of act ive galactic nucleus (AGN) 
activity (e.g.. iFarrah et aTll200ot ). Which of these two con- 
tributions dominates and is primarily responsible for the 
overall dust heating, is still an open question. 

iKewlev et al.l ((2006) presented a robust classification 
scheme of galaxies based on optical emission line ratios of a 
large sample of galaxies from the Sloan Digital Sky Survey 
(SDSS). This scheme allows to discriminate between star- 
bursts, Seyferts, low-ionization narrow emission-line regions 
(LINE Rs), and composit e starburs ts-AGN types . More re- 
cently, lYuan et al.l (|2010T ) used the IKewlev et alj scheme to 
classify a sample of IR selected galaxies, as a function of 
IR luminosity and merger stage. Their results support an 
evolutionary scenario in which ULIRGs are dominated by 
starburst activity at an early merger stage; at intermedi- 
ate stages, ULIRGs would be powered by a composite of 
starburst- AGN activity; and finally, at later stages, an AGN 
would dominate the emission. 

A key feature of ULIRGs is their large dust content, 
which is heated by a central power source, or sources. Since 
optical obscuration is high, radio observations (i.e., extinc- 
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Table 1. The EVN ULIRG sample. Values in columns 2—6 are those provided bv lSanders et al ] [|2003h . The supernova rates in 
column 7 were ob t ained following the empirical relation between CCSN rate and L [8-1000 fim] for starburst galaxies obtained by 
iMattila fc Mc iklc ( 2001), without discounting AGN contribution, if any. In column 8 we show the flux density range a SN would 
display when having a typical type II and type Iln SN peak luminosity, denoted by L ~ irj 27-28 ergs -1 Hz -1 . 



IRAS name 
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on (J2000) 
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Redshift 
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(7) 


(8) 


07251-0248 


07 27 37.5 


-02 54 55 


344 


0.088 


12.32 


5.6 


7-71 


19297-0406 


19 32 22.1 


-04 00 02 


338 


0.086 


12.37 


6.3 


7-73 


19542+1110 


19 56 35.4 


+ 11 19 03 


257 


0.065 


12.04 


3.0 


13-127 


23365+3604 


23 39 01.7 


+36 21 14 


252 


0.064 


12.13 


3.6 


13-132 



tion free) represent the most direct way to distinguish be- 
tween a starburst and an AGN, via the detection of super- 
novae (SNe), supernova remnants (SNRs) and/or compact 
sources at mas-resolution with a high brightness temper- 
ature (Tb), possibly accompanied by a core-jet morphol- 
ogy and usually associated with a high X-ray luminosity. 
Very long baseline interferometry (VLBI) observations have 
been particularly useful, for instance, to discover a popula- 
tion of bright radio SNe and SNRs in the n uclear regions 
of the ULIRGs Arp 22 jSmith et all Il998l ) and Mrk273 
i|Carilli fc Tavlorl |2000| ). This has also been the case for 
LIRGs (Lir > 10 Lq), such as Arp 299 where a prolific 
starb urst and a low-luminosity AGN (LLAGN ) were discov- 
ered (|Neff et al.ll2004l ; |Perez-Torres et alj|20ld . respectively) 
through VLBI observations, or the recent detection of AGN 
activity in a number of LIRGs from the Compact Ob jects 
in Low-power AGN (COLA) sample l|Parra et alj|20ld ). 



2 THE EVN ULIRG SAMPLE 

This is the first of a series of papers presenting European 
VLBI Network (EVN) high-resolution, high-sensitivity im- 
ages of a sample of four of the farthest and brightest ultra 
luminous infrared galaxies (ULIRGs) in the local Universe 
(z < 0.1), part of the project entitled "The dominant heat- 
ing mechanism in the central regions of ULIRGs" (PI: Perez- 
Torres) . 

The sample of ULIRGs we present results from the 
following selection process. We have first selected those 
sources from the IR AS Revised Bright Galaxy Sample 
l|Sanders et alj [20031 ) having logio(im/LQ) > 12, from 
which large supernova rates (u SN ) were expected. We fur- 
ther constrained our sample by selecting those objects with 
5 > — 5° (in order to obtain a good Mv-coverage with the 
EVN), which also appear in the 1.4 GHz Atlas Catalogue o f 
the IRAS Bright Galaxy Sample (|Condon et al.lll990l . [l99r3) , 
as to ensure their radio emission detection. Finally, we se- 
lected those ULIRGs for which neither Multi-Element Radio 
Linked Interferometry Network (MERLIN) nor deep VLBI 
data existed in the literature, and for which MERLIN or 
EVN archival data are not available. The resulting sample 
contains four of the brightest and farthest ULIRGs in the 
local Universe (Table [T]), for which we aimed to unveil their 
dominant heating mechanisms. 

The needed rms to obtain 3a detections of typical 
type II core-collapse SNe (CCSNe) in the most distant 
ULIRGs of our sample, is quite low (~ 2 fiJy, for peak lu- 



minosities ~ 10 27 ergs -1 Hz -1 ) and therefore are well be- 
low our detection limit. Nevertheless, it is also expected 
that more luminous systems provide denser environments, 
which in turn favour the production of very luminous CC - 
SNe (e.g., type Iln SNe in Arp220; iParra et all 120071 ). 
Moreover, the more luminous a radio SN (RSN) is, the 
longer it wi ll take for it to reac h its peak brightness (see 
figure 5 in I Alberdi et al] I2006J ). For instance, the RSN 
AO discovered in th e nuclear region of Arp 299 in 2003 
(see iNeff et al.l [20041 ) . is still det ected after several years 
and remains particularly strong (P erez-Torres et al.l 12009] . 
2010). A similar scenario in the dense nuclear regions 
of our sample of ULIRGs can be expected. Furthermore, 
in the circumnuclear regions of LIRGs and ULIRGs, we 
also expect SN activity to occur. Two remarkable exam- 
ples are SN 2000ft and SN 2004ip, discovered at 600 pc and 
500 p c from the nucl e us of g alaxies NGC 7469 (IColina, et al.l 
l200ll ; lAlberdi et all 120061; IPerez-Torres et al.l l2009h and 



IRAS 18293-3413 (jMattila et all 120071 ; iPerez-Torres etafl 
120071 ). respectively. 

It is worth noting that the empirical relation between 
C CSN rate and L [8-1000 fim \ for starburst galaxies obtained 
bv lMattila fc Meikld <|200ll ): 



yr -1 



2.7 x 10" 12 x 



Lir 
Lq" 



(1) 



assumes no AGN contribution to the IR luminosity. The 
same is true with a similar relation t hat resul t s from the 
combination of equations 20 and 26 in ICondonl (|l992), and 
which yields slightly larger values, i.e., 



y r_1 



3.7 x 10" 



J 



with Lfir = £[40-400 fim] Lq. Thus, if an AGN is present, 
the values for ^ SN in Table [1] represent upper limits, and 
a quantitative estimate of the AGN contribution to the IR 
luminosity is needed before deriving reliable CCSN rates. 



2.1 The case of IRAS 23365+3604 

IRAS 23365+3604 (hereafter IRAS 23365) is t hought to be 
in an advanced merger state (|Sopp et al.lll990T ). There is no 
companion galaxy so far detected with either Very Large 
Array (VLA) or Two Micron All Sky survey (2MASS ) 
observations (see e .g. Sopp et all Il990l : I Yuan et "all 120101 ) . 



iKlaas fc Elsaesserl (|l99ll ) report a companion candidate, a 
small galaxy located ~ 50kpc (projected distance) away 
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from IRAS 23365, which is however not considered to be 
the cause of the extremely irregular and disturbed morphol- 
ogy of IRAS 23365. The optical spectrum of this ULIRG 
seems to be the result of the superposition of LINER and 
Hll-region like components. Such AGN-starburst compos- 
ite spectrum have been confirmed in other studies (e.g., 
IVeron et al1ll997l ; lYuan et al.ll2010T ). Chandra X-ray obser- 
vations have evidenced the pr esence of an AGN ( possibly 
Compton-thick) in this source (jlwasawa et alJfeoill ). 

At a distance of 252 Mpc (1 mas w 1.2 pc), the high 
luminosity of IRAS 23365 (logi O (Li R /L ) = 12.13) corre- 
sponds to a CCSN rate of « 3.6 yr -1 , according to Equation 
[T] As indicated in section this es timate does not con sider 
an AGN contribution. According to lFarrah et al.l |2003r i. the 
AGN contribution in IRAS 23365 is approximately 35 per 
cent of the total Lir, and the rest is due to a starburst, 
from which we infer that v SK w 2.4 yr . 



3 EVN OBSERVATIONS AND DATA 
REDUCTION 

We performed EVN observations of IRAS 23365 quasi- 
simultaneously at L- (y ~ 1.7 GHz or A ~ 18 cm) and C- 
band (y ~ 5 GHz or A ~ 6 cm) in three epochs with a time 
span among them of approximately one year (see Table 

All the epochs were VLBI phase-referenced experiments 
using a data recording rate of 1024 Mbps with two-bit sam- 
pling, for a total bandwidth of 128 MHz. The telescope sys- 
tems recorded both right-hand circular polarization (RCP) 
and left-hand circular polarization (LCP). The data were 
correlated at the EVN MklV Data Processor at JIVE using 
an averaging time of 2 s in the first two epochs, and 4 s in the 
third one, since there was no need for a field of view (FOV) as 
large as 20arcsec. The sources 2134+004 and 3C454.3 were 
used as fringe finders in all the observations. Each epoch 
lasted 6hr, from which a total of ~ 3.7 hr were spent on 
target. Target source scans of 3.5 min were alternated with 
1.5 min scans of the phase reference source. 

The correlated data of every epoch were analysed using 
the NRAO Astronomical Image Processing System (AIPS). 
The overall quality of the visibilities was good and the EVN 
pipeline products were useful for the initial steps of the 
data reduction. To improve the calibration, we edited the 
data to remove artifacts due to radio interference (RFI) 
and included ionospheric corrections where needed. We ex- 
ported the data of all the calibra tors to the Caltech program 
DIFMAP ( Sh epherd et al . 1995) and made images and vis- 
ibility plots of each source. This allowed us to test the per- 
formance of each antenna and to determine gain corrections 
for each. When the gain correction was larger than 10 per 
cent for a given antenna during the whole observing run, we 
applied it to the m>-data using the AIPS task CLCOR. 

In the first two epochs we used J2333+3901 (at 2.9° 
angular distance of target) as a phase reference source. This 
source has a complex structure (see Figure [TJ and it varied 
in flux density at L-band between epochs (see column 7 in 
Table |2j. The subtraction of the phase contribution due to 
the structure of J2333+3901 from the fringe solutions (delay 
and rate) was thus necessary. In spite of this correction, 
the phase referencing of the target source resulted in noisy 
phases. 
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Figure 1. EVN contour maps of the phase reference source 
J2333+3901 at L-band, first | (a) [ and second epoch [(b)] 



In the third epoch we used J2330+3348 (at 3.1° angular 
distance of our target) as a phase calibrator, which being a 
predominantly compact source (at mas angular scales) pro- 
vided a reliable phase reference and calibration. To correct 
the reference position in our first two epochs, and to align 
the three observing epochs, we obtained the shifts in right 
ascension and declination for the first two epochs that make 
their 15<r emission coincide positionally with the 15cr emis- 
sion of the third epoch. We did this by means of the task 
UVSUB in AIPS, in which the data was divided by a point 
source model of 1 Jy at the wanted reference position. 



3.1 Imaging process 

The extended emission of IRAS 23365 is not completely re- 
solved with the available EVN array. The shortest baselines, 
such as Ef-Wb, can recover some of the extended emission. 
If no other short baselines are present (e.g. combinations 
of Jb, Cm and Kn), it is not possible to determine closure 
phases, and the presence of strong sidelobes (of the order of 
the peak) in the dirty map is thus favoured. This situation 
made it very difficult to obtain a reliable image of the target 
source (see the preliminary ma ps of IRAS 23365 presented in 
iRomero-Canizales et ajj 12008). In principle, removing such 
short baselines would solve the problem, at the expense of 
significantly degrading the final image sensitivity. 

To properly map the extended emission, a good cov- 
erage of short baselines (resulting from combinations of at 
least three antennas to determine closure phases) is needed. 
To overcome the lack of short iti>-spacings, a combination 
of Gaussian model fitting and imaging algorithms can be 
used. This is a widely used method for mapping the struc- 
ture of outflows at VLBI scales (see e.g.. iRastorgueva et al.l 
2011), specially for the cases in which faint diffuse emission 
is present together with the bright compact one. Epochs 
1 and 2 were affected by poor short-baseline Mii-coverage 
(the Cm-Kn baseline had a severe amplitude problem and 
it was not used). In epoch 3, baselines Jb-Kn, Jb-Cm, Cm- 
Kn and Ef-Wb were present, thus permitting to determine 
closure phases for the short baselines. As a result, no strong 
sidelobes affected the imaging process at this epoch. Nev- 
ertheless, for this epoch we also used a Gaussian model fit- 
ting combined with clean components in order to obtain 
consistent results with those of the first two epochs. This 
was done within the C altech imaging programme DIFMAP 
l|Shepherd et al.lll995r i. We exported the resulting images 
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Table 2. Parameters of the EVN observations. The stations (location, diameter) used in the different observing runs are: Ef- 
Effelsberg (Germany, 100m), Wb-Westerbork array (NL, 14x25m), Jbl-Lovell (UK, 76m), Jb2-MKII (UK, 25m), On-Onsala 
(Sweden, 25 m), Mc-Medicina (Italy, 32 m), Nt-Noto (Italy, 32 m), Tr-Torun (Poland, 32 m), Ur-Urumqi (China, 25 m), Cm- 
Cambridge (UK, 32m), Kn-Knockin (UK, 25m), Ar-Arecibo (Puerto Rico, 305m), Ys-Yebes (Spain, 40m). In column 6 we show 
the phase reference source used in each epoch and in column 7 their associated peak intensities. 



Label 


Project 


Observing 


Frequency 


Participating 


Phase 


^peak 






date 


(GHz) 


stations 


calibrator 


( Jy/beam) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


LI 


EP061A 


2008-02-29 


1.7 


Ef, Wb, Jbl, On, Mc, Nt, Tr, Ur, Cm 


J2333+3901 


0.34 ± 0.02 


CI 


EP061C 


2008-03-11 


5.0 


Ef, Wb, Jbl, On, Mc, Nt, Tr, Ur, Cm 


J2333+3901 


0.21 ± 0.01 


L2 


EP064D 


2009-03-07 


1.7 


Ef, Wb, Jb2, On, Mc, Nt, Tr, Ur, Cm, Kn, Ar 


J2333+3901 


0.48 ± 0.02 


C2 


EP064B 


2009-02-28 


5.0 


Ef, Wb, Jb2, On, Mc, Nt, Tr, Ur, Cm, Kn, Ar 


J2333+3901 


0.23 ± 0.01 


L3 


EP064J 


2010-03-08 


1.7 


Ef, Wb, Jbl, On, Mc, Nt, Tr, Ur, Cm, Kn 


J2330+3348 


0.65 ± 0.03 


C3 


EP064L 


2010-03-20 


5.0 


Ef, Wb, Jbl, On, Mc, Nt, Tr, Ur, Cm, Kn, Ys 


J2330+3348 


0.64 ± 0.03 



Table 3. EVN observational data, as measured from the images shown in Figure[2] Columns (2) and (3) Peak position coordinates 
given with respect to a( J2000) = 23 h 39 m 01 s 0000 and <5(J2000) = 36°21'08'.'000. The errors in position (within parentheses), given 
in mas, were estimated by adding in quadrature the random errors in the position of the target (tar) and the phase reference 
calibrator (ref): v /(FWHM tar /(2 X SNRt^)) 2 + (FWHM rcf /(2 X SNR rci )) 2 , where SNR is the signal to noise ratio, and FWHM 
was taken as the projection of the beam major axis on both a and <5 axes. Column (4) rms noise in the maps. Column (5) Peak 
intensities. Column (6) Flux densities measured in regions enclosing 5cr level of the emission. Column (7) Matched C-band flux 
densities, covering the 5<r level L-band emission region; this is done by adding the C-band Stot to a 3<r emission in the remaining 
part. The uncertainties for the measurements shown in columns 5-7 have been estimated by adding in quadrature the rms noise in 
the map plus a 5 per cent uncertainty in the point source calibration. Column (8) Characteristic size (in both a and <5) occupied 
by the emission at a 5cr level.) 



Label 


Aa(J2000) 


A<5(J2000) 


rms 


•Speak 


Stot 


^match 


R a x R§ 




( s ) 


(") 


( /xjy/beam) 


( /^Jy/beam) 


(mjy) 


(mjy) 


(pc 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


LI 


0.2600 (0.5) 


0.592 (0.5) 


28 


786 ± 48 


7.99 ± 0.40 




207 x 221 


CI 


0.2614 (0.7) 


0.603 (0.7) 


16 


303 ± 22 


0.32 ± 0.02 


1.33 ± 0.07 


68 x 69 


L2 


0.2616 (0.7) 


0.598 (0.7) 


25 


466 ± 34 


5.42 ± 0.27 




221 x 220 


C2 


0.2607 (0.5) 


0.603 (0.5) 


23 


584 ± 37 


1.11 ± 0.06 


2.38 ± 0.12 


98 x 94 


L3 


0.2615 (0.6) 


0.566 (0.6) 


30 


640 ± 44 


8.54 ± 0.43 




241 x 259 


C3 


0.2608 (0.3) 


0.599 (0.3) 


18 


875 ± 47 


2.45 ± 0.13 


3.52 ± 0.18 


111 x 127 



back into AIPS to analyse them and to produce the final 
maps that we present here (see Figure [2]) . 

In Table [2] we show the stations that participated in 
each observation. For different reasons, we lost some anten- 
nas and/or baselines and the final images were produced us- 
ing the visibilities resulting from slightly different arrays. For 
instance, in the second epoch we lost Ur and Ar, and thus 
the resolution was compromised by the loss of the longest 
baselines. On the other hand, in the third epoch we had a 
good coverage of the short baselines (from combinations of 
Cm, Kn, Jb, Ef and Wb), which eased the reconstruction 
of the extended emission. To allow comparisons among the 
different epochs and frequencies, we used the same convolv- 
ing beam (that from the epoch with the worst resolution: 
26x38 mas 2 at 46°) for the imaging process and sampled 
the beam using the same cell size (= 4 mas), and natural 
weighting for all epochs. The resulting images for the three 
EVN epochs at the two different frequencies (1.7 and 5 GHz) 
are shown in Figure [2] The actual array used in the different 
epochs, is shown in a label at the upper right corner of each 
image. In Figure [3] (bottom) we also show the third EVN 
epoch at both 1.7 and 5 GHz, as imaged with the natural 
beam of the observation at 1.7GHz (19x25 mas 2 at 30°). 



4 MERLIN AND VLA OBSERVATIONS 

Simultaneously with our second EVN epoch, we also ob- 
served IRAS 23365 at both L- and C-bands with MERLIN 
(see Table [4|, including the following antennas: Defford, 
Cambridge, Knockin, Darnhall, Mark 2 and Pickmere, ob- 
serving with a bandwidth of 15 MHz (in both circular polar- 
isations). OQ208 was used as amplitude calibrator (1.1 Jy 
at L-band and 2.5 Jy at C-band) and J2333+3901 (0.8 Jy 
at L-band and 0.34 Jy at C-band) as phase calibrator. For 
phase-referencing, duty times of 7min/1.5min in L-band, 
and 2.5 min/1.5 min in C-band were used, for a total time 
on source of 4 and 2.5 hr at each band, respectively. 

We also analysed archival VLA (A-configuration) 
data at L- and C-ba nds (project: AB660, reported in 
Baan & Klocknerl lioblSI) to compare with the MERLIN and 
EVN images. The observations were performed with a band- 
width of 50 MHz (in both circular polarisations). 3C48 
(16.0 Jy at L-band and 5.4 Jy at C-band) was the flux cali- 
brator and 0025+393 (0.7 Jy and 0.6 Jy at L- and C-band, 
respectively) the phase calibrator, which is found at 9.6° 
angular distance of the target source. 

We followed standard procedures within AIPS for the 
data reduction. Details on the VLA and MERLIN observa- 
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tions are shown in Tableland the resulting images in Figure 
[3] We used matched baselines (in wavelengths) to obtain the 
images for the two different frequencies of each array, to en- 
able the comparison of information at the same scales. For 
the VLA images we used a common uif-range of 11.5 to 162.7 
kA and same convolving beam of 1.3 x 1.3arcsec 2 , while we 
restricted MERLIN images to 112.6 to 1191.8 kA, that re- 
sulted in a convolving beam of 0.10x0.25 arcsec 2 , at 30°. We 
did not perform any tti>-restriction in the case of the EVN 
data, to optimize the sensitivity and ui>-coverage for each 
observing epoch. 



5 RESULTS 

IRAS 23365 has been observed at different resolutions (EVN, 
MERLIN and VLA) and at two different frequencies. 
This allows a comparison among the different linear scales 
mapped with different arrays. In the following, we present 
our results regarding morphology, radio emission, radio spec- 
trum and magnetic field of IRAS 23365. The different param- 
eters measured from the three epochs of EVN observations 
(see Figure [2} are presented in Table and measurements 
from the VLA archival data and the MERLIN observations 
are presented in Table [3] In Table [5] we show the estimates 
from the measurements at different scales as shown in Figure 



m 



5.1 The IRAS 23365 structure: from kpc- down to 
pc-scales 

The radio images of IRAS 23365 shown in Figure [3] cover 
the structure of this galaxy at three different scales: galactic 
(with the VLA), circumnuclear (with MERLIN) and nuclear 
(with the EVN). We use the third epoch of EVN observa- 
tions to compare with the VLA and MERLIN images, since 
that was the epoch with the best compromise between an- 
gular resolution and sensitivity (see Section T3. II for details). 
For doing this, we re-imaged the third EVN epoch (L- and 
C-band maps) using as convolving beam that one which re- 
sulted from the L-band (19x25 mas 2 at 30°, see bottom of 
Figure O. 

At galactic scales, the emission at both L- and C- 
bands is unresolved and appears concentrated in a zone of 
~ 4.5 kpc in size. At circumnuclear scales, the emission is 
concentrated in the inner 0.5 kpc region and displays some 
extended structure on top of an unresolved component. 

At the highest resolution in L-band, the nuclear region 
has a size >200pc at all the EVN epochs (see Table [3]), 
and shows variations in its morphology (see Figure [2]). At 
C-band, the emitting region is about 100 pc and its structure 
remains quite compact in the first two epochs, whilst some 
more extended emission is traced in the third epoch (see 
Figure [2J • A single Gaussian fit is inaccurate for obtaining 
the deconvolved size of the emitting region, at least for the 
emission at L-band, due to the wealth of extended emission. 
We thus characterize the area of the emitting region with 
the size of the source in both the right ascension and dec- 
lination axes, R a and Rs, respectively (Table [3]). There are 
some features outside the nuclear regions that, while having 
peaks slightly above 5<r in our second epoch (Figures |2(d)| 
and 2(e) I, are not seen neither in our first epoch, nor in our 



third observing epoch. While these could be real features (in 
particular, the compact source detected at both frequencies 
with Act, AS ~ 150 mas), we conservatively consider them 
as tentative detections (see Section [3~T1 for details) and there- 
fore are not discussed here. 

We note that the size of the emission area increases with 
time through the different EVN epochs at both frequencies 
(see Table [3]), and also displays different morphology (see 
Figure [2]), especially at L-band. Whereas sensitivity does 
not seem to vary drastically among epochs, the observa- 
tions were performed at different hour angles and thus the 
Mu-plane was sampled at different orientations. Hence, the 
differences in size and morphology could have been affected 
by the different u^coverages. 

Regardless of the used array (i.e., VLA, MERLIN or 
EVN), and albeit of using matched baselines (at least for 
VLA and MERLIN), the emission in L-band consistently 
occupies a larger extension than that at C-band, around a 
factor of 2 in the case of the EVN, as seen in Figure[3] where 
we show for comparison the VLA, MERLIN and EVN (third 
epoch) images. This can be explained by the longer lifetime 
of accelerated electrons emitting synchrotron radiation at 
lower frequencies (see Section [5. 4p . 

We also note that the peak positions at the two dif- 
ferent frequencies are not coincident neither for MERLIN 
nor for the EVN. In the case of the VLA, we do not have 
the required angular resolution to confirm any shift; how- 
ever, at the higher resolution provided by both MERLIN 
and EVN, a shift of the C-band peak towards the North- 
East direction is evident, while that at L-band is shifted 
towards the South- West (see Tables [3] and [4] and Figure [3]). 
This result is consistent for all the epochs and at the dif- 
ferent angular resolutions provided by EVN and MERLIN. 
We interpret those shifts of the emission peaks as evidence 
for at least two different populations of radio emitters being 
present in the nuclear region (see Section \5. 2 [) . Furthermore, 
the peak component is variable both in position and in in- 
tensity among EVN epochs, and in each epoch, being also 
different between frequencies. These facts give evidence of 
the source variability within the innermost nuclear region. 



5.2 The radio emission and radio spectrum at 
different scales 

We mentioned in the previous section that the radio emis- 
sion at different frequencies seen at the different resolutions 
(except perhaps for the VLA), peaks at different positions. 

Thus, a peak spectral index defined as S pea k oc v pcak would 
be meaningless. We use instead the peak of the pixel-to-pixel 
spectral index distribution (a") as obtained with AIPS. 

The radio emission of both galactic and circumnuclear 
regions of IRAS 23365 mapped with the VLA and MER- 
LIN, respectively, is stronger at L-band than at C-band (see 
columns 7 and 8 in Table [4}. Consequently both total spec- 

L-C 

tral indices (Stot oc v * ot ) and peak pixel-to-pixel spectral 
indices (OpI°) are steep, as shown in columns 4 and 5 of Ta- 
ble[5] Steep spectral indices are an indication of non-thermal 
emission. We note however that the Tb values (column 2 in 
TableJS]) at galactic (L- and C-bands) and circumnuclear (C- 
band) scales, are in principle consistent with either thermal 
emission, or with synchrotron emission suppressed by free- 
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Table 4. Parameters and observational data of complementary VLA and MERLIN observations (see the corresponding contour 
maps in Figure[3]l. The observations are labelled with a two- letter code, in which the first letter corresponds to the array (V=VLA 
or M=MERLIN), and the second letter to the observed frequency band (Column 3). Columns (4) and (5) Coordinates of the peak 
position, given with respect to a(J2000) = 23 h 39 m 01?0000 and 5(J2000) = 36°21'08'.'000. Column (6) - rms noise in the maps. 
Column (7) Peak intensities. Column (8) Flux densities measured in regions enclosing 5it level of the emission. The uncertainties 
for the measurements shown in columns 4-5 and 7—8 have been estimated as described in Table [3] Column (9) Major and minor 
axes, obtained by fitting a Gaussian to the source. Column (10) Characteristic size (in both a and <5) occupied by the emission at 
a 5<r level. 



Label 


Observing 


Freq. 


Aa(J2000) 


A<5(J2000) 


rms 


'Speak 


Stat 


0M X m 


R a x R$ 




date 


(GHz) 


H 


(") 


( /ijy/beam) 


( mjy/beam) 


(mjy) 


(arcsec 2 ) 


(kpc 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


VL 


1992-12-14 


1.4 


0.252 (6.1) 


0.54 (6.1) 


180 


19.14 ± 0.97 


25.18 ± 1.27 


0.93 x 0.58 


4.48 x 4.38 


VC 


1992-12-14 


4.9 


0.261 (3.3) 


0.59 (3.3) 


50 


9.97 ± 0.21 


10.82 ± 0.54 


0.46 x 0.30 


4.54 x 4.16 


ML 


2009-03-06 


1.6 


0.260 (1.6) 


0.55 (4.0) 


200 


6.29 ± 0.37 


13.90 ± 0.72 


0.17 x 0.14 


0.36 x 0.66 


MC 


2009-02-25 


5.0 


0.264 (3.6) 


0.56 (8.9) 


170 


2.39 ± 0.50 


5.20 ± 0.31 


0.17 x 0.16 


0.32 x 0.54 



Table 5. Physical quantities of IRAS 23365 estimated from VLA, MERLIN and EVN observations. We have used the measured values of 
Speaki ^tot and the largest linear size between R a and Rg, from Tables [3] and [4] In the case of the EVN, we only include the third epoch of 
observations, re-imaged with a convolving beam of 19x25 mas 2 at 30° (see Figure[3j|. We consider the emission at both bands as measured 
from the region delimited by the 5<r L-band emission, i.e., using Stot for L-band, and S ma t c h for C-band, and the same linear size for both. 
Column (2) Brightness temperature calculated using as the source' solid ang le Q s = 7r(41og2)- 1 (6 M x 6 m ) for the VLA and MERLIN 
(Table|4]l; for the EVN no reliable Gaussian fit can be made to the source (especially in the L-band map) in order to find its deconvolved size, 
and we have used instead f2 s = 7r(41og2)~ 1 (FWHM M X FWHM m ), where FWHM M and FWHM m are the major and minor synthesized 
beam fitted FWHM. Column (3) Monochromatic luminosity at the frequency v. Column (4) Two-point peak pixel-to-pixel spectral index 
determined from the pixel-to-pixel spectral index distribution obtained with AIPS. Column (5) Two-point total spectral index determined 
between the L- and C-band flux densities (Stot oc u a t°t ). Column (6) Integrated isotropic radio luminosity, considering that the spectral 
behaviour of the source is straight between L- and C-band. Column (7) Equipartition magnetic field. Column (8) Characteristic lifetime of 
electrons subject to B cq , undergoing radiative synchrotron losses. 



Label logT B L v o£;° ofc£ Lr B cq r syn 

(K) (10 29 ergs" 1 Hz" 1 ) (10 39 ergs" 1 ) (/|G) (Myr) 

(1) (2) (3) (4) (5) (6) (7) (8) 



VL 


4.45 ± 0.02 


19.11 ± 0.97 


VC 


3.61 ± 0.02 


8.21 ± 0.41 


ML 


5.42 ± 0.02 


10.55 ± 0.55 


MC 


3.96 ± 0.03 


3.95 ± 0.24 


EL 


6.81 ± 0.02 


4.94 ± 0.25 


EC 


5.57 ± 0.02 


2.58 ± 0.13 



-0.46 ± 0.06 -0.69 ± 0.06 4.1 
-0.66 ± 0.13 -0.89 ± 0.07 2.0 
0.52 ± 0.13 -0.59 ± 0.06 1.1 



18.4 
18.1 
77.1 
90.8 
174.7 
174.7 



11.4 
6.2 
1.1 
0.6 
0.4 
0.2 



free absorption from e.g., HII regions. The calculated value 
for the free- free opacity (rg) implies that thermal emission 
should be optically thick. Therefore, the bulk of emission at 
L- and C-bands corresponds to optically thin non-thermal 
synchrotron emission, slightly affected by free-free absorp- 
tion at galactic (rg ~ 0.02) and circumnuclear (rg ~ 0.20) 
scales. 

Regarding the nuclear region (mapped with the EVN), 
the large Tb values are consistent with pure non-thermal 
emission. The attained angular resolution and the presence 
of strong extended (a few mjy; see Tables [3l and Q radio 
emission, prevents us from directly detecting individual faint 
(see Table [TJ compact sources, e.g. SNe. However, we note 
that .Speak and Stot show variations at both frequencies dur- 
ing our EVN monitoring campaign (see Table [3} . Whereas 
in C-band, 5 pC ak and 5*tot increase with time, in L-band 
these diminish in the second epoch, and then increase in the 
third one, thus indicating the variability of sources and/or 
the appearance of new ones within the nucleus, e.g., new 
SNe, accounting to the expected SN rate (« 2.4yr _1 ). This 
non-correlated behaviour at both frequencies is indicative of 



nuclear activity that becomes transparent first at C-band 
and later at L-band. 



5.3 Spectral index distribution at mas-scales 

Let us now concentrate in the spectral indices correspond- 
ing to the EVN images. Considering the total flux densities 
as measured from the region within the 5a L-band emis- 
sion (i.e. L-band flux from column 6, and C-band flux from 
column 7 in Table [3} , the total spectral indices (a^ f) are 
steep for all the epochs. However, the situation is different 
for the peak pixel-to-pixel spectral index (apTp), which is 
evolving with time. The distribution of apTp is shown in 
Figure [4] for the three EVN epochs. In the first epoch of 
EVN observations, apTp is steep, then it becomes inverted 
in our second epoch, and starts to decrease (although be- 
ing still inverted) in the third epoch to presumably become 
steep again. This is clear evidence of the variation in flux 
of sources within the innermost nuclear regions, and/or ap- 
pearance of new sources (e.g., SNe) which would be seen 
first at higher frequencies and later on at lower frequencies 
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|Weiler et alj|2002h . in agreement with our results. We also 
note that for the three EVN epochs, apTp (which is given 
pixel by pixel as shown in Figure [4]), becomes steeper as 
measured towards the edges of the C-band emission, where 
the noise at C-band starts to dominate, whilst there is still 
extended emission detected at L-band. This is a consequence 
of the ageing of the population of electrons radiating syn- 
chrotron emission (see Section f5.4|) . 



5.4 The magnetic field in the energy budget of 
IRAS 23365 

In previous sections, we have gathered information about 
the ongoing non-thermal activity of IRAS 23365 at differ- 
ent scales. In a ULIRG environment, we expect SNe, SNRs 
and/or an AGN to be the engines responsible for producing 
high energy particles which will interact with the galactic 
magnetic field, thus ge nerating synch rotron radiation (dom- 
inating at v < 30 GHz; ICondon|[l992l ). The energy thus pro- 
duced, will be present in the form of relativistic particles 
and magnetic field. In the following, we investigate the en- 
ergy budget (due to synchrotron radiation) of IRAS 23365 
at different scales, i.e., as estimated from observations with 
different arrays (EVN, MERLIN and VLA). We only con- 
sider the third epoch of observations with the EVN, to com- 
pare with the results from the VLA and MERLIN, since this 
epoch was the one which had the least imaging problems (see 
Section 13. ip . 

We can estimate the average equipartition magnetic 
field based on the ra dio emission of IRAS 23365 as follows 
(sec |Pacholczvklll970h . 



4> 



(lO?) (lkpc) ( 



10 39 



(2) 



where <j> is the filling factor of fields and particles, k is the 
ratio of heavy particle energy to electron energy, and c 12 is a 
function that depends on the minimum and maximum fre- 
quencies considered, and of the two-point spectral index, 
a tot , which is estima ted based on those two frequencies 
(see I Pacholcz vkl ll970T l . Lr is the integrated isotropic radio 
luminosity between the minimum and maximum frequen- 
cies used, and R is the linear size occupied by the emis- 
sion, taken as the larger value between R a and Rs in each 
case (see column 10 in Table [4] for the VLA and MER- 
LIN). For the third EVN epoch, we determined a maxi- 
mum linear size R « 0.2 kpc with TVDIST within AIPS. 
For simplicity, we consider <j> = 0.5 and k — 100 (see e.g., 
iPerez- Torres fc Alberdill2007l ). 

In Table [5] we show the average values for al a t , Lr and 
Boq, obtained within the emission regions sampled by the 
different instruments. In the innermost nuclear region (im- 
aged with the EVN), the strength of the magnetic field is 
larger than the one measured at lower resolutions. This is 
expected, since the plasma in the central regions should be 
denser than in the outer regions, and thus the magnetic field 
lines therein, frozen within the plasma, should be more con- 
centrated. The average magnetic field under energy equipar- 
tition, for the emission measured with the VLA, MERLIN 
and the EVN, would be 18, 84 and 175 ^G respectively. The 
latter value represents the peak of L? cq coming from the very 
central region. If the synchrotron spectrum holds beyond the 



C-band frequencies, e.g. to 20 GHz, the estimated values for 
S eq , would only be ~ 15 per cent larger. 

Our obtained L? cq value at galactic scales is consistent 
with that of a galaxy in advanced interaction state, prob- 
ably close to nuclear coalesce nce, according t o VLA stud- 
ies of interacting galaxies by iDrzazga et al.l (|20 1 lh . Like- 
wise, the Beq value at nuclear scales is similar to that 
found th rough VLBI studies of the ULIRG IRAS 17208-0014 
(144 pG; iMomiian et alj 2003), which is also an advanced 
merger. 

Consid e ring t he obtained values for -Beq, and following 
|Pacholczvkl(|l970h , we can calculate the lifetime of the elec- 
trons with energy L m i n , which move in a magnetic field of 
strength L? eq , thus emitting synchrotron radiation around a 
critical frequency v c . This is, 



'syn — R2 , 

(where c 2 is a constant) 



Myr 



1.06 x 10 d 



-, -1/2 



(3) 



From Section 15.11 we know that the radio emission at 
L- and C-bands has a different extent and peaks at different 
positions within the nuclear region, which strongly suggests 
the presence of different populations of particles. This is 
more evident in the nuclear region mapped with the EVN: 
in the innermost region, where there is an overlap between 
the emission at the two different frequencies, there would 
be a concentration of very energetic, short-lived particles, 
whereas the outer region would be populated by less ener- 
getic, long-lived particles, which had had time to diffuse from 
the inner regions into the outer ones. In Table[5]we show the 
values for r syn , assuming that the critical frequency is either 
that of the L-band or the C-band. In all cases we see that the 
L-band emission is tracing the emission from an older popu- 
lation of electrons (regardless of the resolution) than the one 
emitting at C-band frequencies. The putative AGN together 
with an ensemble of SNe, for which evidence has been found 
in other studies (see Section 12. ip , must be located within 
the C-band emission region as seen with the EVN, where 
the magnetic field strength is larger, and where a composite 
spectrum (which varies with time) has been found (Section 

El. 

We note that the radio lifetime of the emitting source 
is not only determined by r syn . The lifetime of relativistic 
electrons might also be affected by Compton losses given by 

25.2 

TC 



ffollowing |pacholczvklll970r ) since the electrons are immersed 
in a radiation field, 

47T L b ol 



U, 



for which we take Lnt as a good approximation to the 
bolometric luminosity Lboi- U va d varies from ~ 2.5 x 
10 -7 ergcm~ 3 at galactic scales, up to ~ 2.8 x 10 -4 erg cm -3 
at nuclear scales. To compare the different losses, we calcu- 
late their ratio, 



t c /t 



25.202^ 
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Table 6. Parameters derived for the compact components ECl, ELI and EL2 found within the extended emission in the third 
epoch of EVN observations (Figure[3]l. Columns (2) and (3) Coordinates of the peak position, given with respect to a(J2000) = 
23 h 39 m 01?0000 and 5(J2000) = 36°21'08'.'000. Column (4) Peak intensity corrected by the zero level emission. Column (5) 
Monochromatic luminosity. Column (6) Brightness temperature, considering that the solid angle subtended by the source is that 
subtended by the synthesized beam (as in Table [5} . 



Label 
(1) 


Aa(J2000) 

( s ) 
(2) 


A<5(J2000) 

(") 
(3) 


^pcak 

( /xjy/beam) 
(4) 


L v 

(10 28 ergs" 1 Hz" 1 ) 
(5) 


logT B 
(K) 
(6) 


ECl 
ELI 
EL2 


0.2608 (0.4) 
0.2597 (1.7) 
0.2616 (1.4) 


0.599 (0.4) 
0.586 (1.7) 
0.561 (1.4) 


352 ± 35 
150 ± 32 
184 ± 33 


2.67 ± 0.26 
1.14 ± 0.24 
1.39 ± 0.25 


4.56 ± 0.04 
5.15 ± 0.09 
5.23 ± 0.08 



We note that at all scales (nuclear, circumncuclear and 
galactic) and at both L- and C-band, we obtain r syn 3> tc, 
with a ratio ranging between 6.5 x 10 -6 (nuclear scales) and 
8.8 x 10~ 5 (galactic scales); i.e., the energy density of the 
radiation field, greatly exceeds the magnetic energy density. 
Radiative losses are important and we argue that there is 
need for injection of new electrons or a continuous acceler- 
ation to halt the energy depletion, otherwise radio emission 
would not be visible. 

The re-acceleration or injection of new electrons in a 
(U)LIRG environment, is very likely provided in SN-shells 
and SNRs by first order Fermi acceleration. The presence 
of SNe, SNRs and a strong magnetic field in IRAS 23365, 
agrees with this scenario. 



5.5 The nuclear region in the third EVN epoch 

Among the EVN observing epochs, the third one benefited 
from a better «?>coverage, and thus resulted in a smaller nat- 
ural beam (19x25 mas 2 at 30° at L-band). The L-band map 
(Figure [3] bottom-left) shows the presence of two compact 
sources (ELI and EL2) within the nuclear region, without 
counterparts at C-band. On the other hand, the compact 
source that dominates the emission at C-band, labelled as 
ECl (Figure [3] bottom-right), has no compact counterpart 
at L-band, although extended emission is present. 

To obtain ECl, ELI and EL2 peak intensities, we 
first estimated the background emission where these com- 
pact sources lay. We solved for the 'zero level' emission 
(Sb g ~ 179 ± 24 ^Jy at C-band and S bg ~ 211 ± 12 ^Jy at 
L-band) using the task IMFIT within AIPS. We then sub- 
tracted this value from the maximum intensity found at the 
positions of each compact source, in order to obtain their 
5 P eak- In Table [6] we give the positions for ECl, ELI and 
EL2, their peak intensities, as well as their estimated L„ 
and Tb, which are indicative of a non-thermal origin. 

ECl lays on a region where apTp changes with time, sug- 
gesting variability within this region. We note that ECl is 
confined to a small area in the first epoch, and then appears 
to increase in size, as we have mentioned in Section T5. II ELI 
lays in a region with basically no C-band emission, unlike 
the region where EL2 lays where more extended emission is 
being traced from the first epoch to the third one at C-band. 
As a consequence, ELI lays in a region which maintains a 
very steep through time, whilst EL2 is found in a re- 
gion with varying apTp (Figure |4}. However, we cannot rule 
out that the apTp variations at EL2 are intrinsic. 



In the case of ECl, both the variability of the radio 
emission (see Section 15. 2p and of the spectral index distri- 
bution (see Section r5.3l and Figure|4]) are indicative of recent 
non-thermal activity (probably due to SNe and/or AGN ac- 
tivity). ELI and EL2 display brightness temperatures simi- 
lar to those expected from either type II SNe or SNRs. We 
note that the maximum linear size for ELI and EL2 is set 
by the beam size to ~ 30 pc, which is too large for char- 
acterising either an individual SN or a SNR. A scenario in 
which ELI and EL2 are clusters of SNe is difficult to recon- 
cile with the absence of peaks of emission at C-band in all 
the EVN observing epochs, and with the behaviour of apTp 
at both ELI and EL2. These facts suggest that there is no 
recent activity from young SNe in those regions, and favour 
an scenario in which ELI and EL2 are dominated by an old 
population of radio emitters. 



6 SUMMARY AND DISCUSSION 

We have presented state-of-the-art radio interferometric 
images of IRAS 23365, one of the brightest and farthest 
ULIRGs in the local Universe (z < 0.1). 

Our images reveal the presence of a nuclear region, pos- 
sibly a starburst-AGN composite, with an approximate size 
of 200 pc in L-band, and about 100 pc in C-band. We find 
that the L- and C-band radio emission peak at different po- 
sitions, thus suggesting that the nuclear region is composed 
of at least two zones, dominated by distinct populations of 
radio emitters. 

In the region where the L- and C-band emission overlap, 
there is evidence for ongoing non-thermal activity, charac- 
terised by very energetic, short-lived particles . During our 
EVN monitoring of IRAS 23365, we have found flux den- 
sity variability in the overlapping region, thus resulting in 
a variation of the spectral index. This can be explained by 
the flux density variations of sources therein (SNe, AGN, 
etc.) and/or by the appearance of new sources (e.g., SNe) 
which would be seen first at higher fre quencies and later 
at lower frequencies (|Weiler e t al. 2002j). The edges of the 
overlapping region characterised by less energetic, long-lived 
particles, would be dominated by an old population of radio 
emitters, probably clumps of SNRs, for which we have found 
two candidates in the third L-band EVN epoch. These facts 
agree with the c lassification of IRA S 23365 as a composite 
system, made bv lYuan et al.l l|2010h . 

The radio source lifetime at different scales (as seen 
with the VLA, MERLIN and the EVN arrays) and at both 
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L- and C-bands, is limited by Compton losses. The SNe and 
SNRs, for which we have found evidence, are likely providing 
the mechanism of re-acceleration, or replenishment of new 
electrons that is needed to halt the radio energy depletion. 

We have found that the equivalent magnetic field 
strength at galactic (mapped with the VLA) and nuclear 
scales (mapped with the EVN), 18 and 175 fiG, respectively, 
correspon d to that of a galaxy in an advanced stage of in- 
teraction ijDrzazea et al.l l201ll : iMomiian et all 12003). The 
magnetic field in both nuclear and circumnuclear regions 
is stronger than at galactic scales, thus implying that the 
lifetime of the electrons undergoing synchrotron losses is 
shorter ( < 1 Myr) in the innermost nuclear regions (with 
linear size^R < 0.5 kpc) of IRAS 23365, and larger ( > 1 Myr) 
in the outer regions (R > 4 kpc) . 

Our study of IRAsli3365 (at z ~ 0.06) has shown that 
high-resolution, high-sensitivity observations are needed if 
we are to make significant improvement in the detailed un- 
derstanding of nuclear and circumnuclear starbursts in the 
local Universe. The resolution we attained using a maximum 
baseline length of approximately 7,000 km, is not enough to 
resolve individual compact sources (e.g., SNe, SNRs, AGN) 
from each other, within the nuclear region of IRAS 23365; 
yet, it could be possible to infer the activity of such com- 
pact sources, by carefully monitoring variations of total 
flux density and spectra l inde x distribution. For instance, 
iRomero-Canizales et al.l (|201 ll ) were able to directly detect 
SN activity in the Bl-nucleus of Arp299, by carefully mon- 
itoring the variations in Stot , over several years of VLA ob- 
servations. Without spatially resolving each individual SNe, 
they estimated a lower limit for the Vsn in that LIRG. In 
the case of IRAS 23365, where a larger number of SNe are 
expected each year, several observations per year would be 
needed to perform such an indirect study of the SN pop- 
ulation in its nuclear region, provided that we are able to 
distinguish between AGN outbursts and SN explosions. 

IRAS 23365 is a good example of the situation to be 
faced when observing galaxies at higher redshifts. It is ex- 
pected that the Square Kilometer Array (SKA), with a max- 
imum baseline length ~ 3, 000 km, will allow the detection 
of source s as faint as 50 nj y, e.g. CCSNe, exploding at 2 ~ 5 
(see e.g.. lLien et alJfeoilT l. However, the angular resolution 
will be a strong limiting factor. In those cases where the nu- 
clear and even the circumnuclear regions (i.e., where we ex- 
pect most of the SN activity to occur) of the host galaxy can- 
not be resolved out into their different components, SKA's 
high sensitivity might be of great use to indirectly detect SN 
activity through the monitoring of flux density variations. 
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Figure 4. Evolution of the spectral index distribution (5 p i xc i ~ 

L-C 

v p-p ), given pixel by pixel (as obtained from AIPS) in the 
nuclear zone, which we expect to be composed by non-thermal 
sources (e.g. SNe and AGN). Note that while the peak of the 
spectral index distribution varies in the very central region, it 
becomes steeper towards the edges, thus representing older radio 
emission. The sources EC1, ELI and EL2 are marked with crosses 
(see Section [531 . 
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Figure 2. IRAS 23365 at L-band (left: |(a)| |(d)| [(g)) , C-band (middle: | (b) [ [(e)] |(h) | , and L-band contours overlaid on grey scale C-band 
images (right : [(c)||(f)[[(i)| |, in three different epochs (top to bottom), with the grey scale in /xjy ranging from 3tr value of the noisiest map 
in C-band (second epoch, i.e., <r = 23 /ijy/beam), up to the peak intensity value of each epoch at C-band. Dashed contours represent 
— 3<t levels. All the images have been degraded to the epoch with lowest resolution (L2, following the labels in Tabled, being therefore 
convolved with the same beam size: 26x38 mas 2 at 46°. The maps are centred at 23 h 39 m 01?29, +36 o 21'08'.'59 (J2000). At the top of 
each image, we provide the EVN project code, observation date, observed wavelength and the list of antennas used for mapping. In the 
lower left corner we show the peak intensity, noise and contour levels information. The size of the nuclear region seems to increase from 
one epoch to another at both frequencies, although images from the same frequency were convolved with the same beam. We consider 
that this is an effect of the difference in the arrays used, rather than being an intrinsic change in the source, albeit this possibility cannot 
be neglected. It is worth noting that the size of the nuclear zone is consistently larger at L-band than its counterpart at C-band. This 
can be explained by the longer lifetime at lower frequencies of the electrons being accelerated in the innermost nuclear regions. 
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Figure 3. IRAS 23365 as seen at different frequencies and resolutions (higher resolution from top to bottom; and higher frequency from 
left to right, see labels). All the maps are centred at 23 h 39 m 01?29, +36°21'08y59 (J2000). Dashed contours represent -3<r levels. The 
cross in the EL image indicates the peak position (ECl) as measured from the EC image, and the crosses in the EC image indicate the 
positions of components ELI and EL2. The non-coincidence of the peaks at the two different frequencies (see also ML and MC images), 
indicates the presence of distinct population of sources in the innermost nuclear region. We note that at C-band, IR AS 23365 remains 
practically unresolved, whilst L-band images show more extended structure. 
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